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ABSTRACT 

Massive Black Holes detected in the centers of many nearby galaxies show an approximately linear 
relation with the luminosity of the host bulge, with the black hole mass being 0.001-0.002 of the bulge 
mass. Previous work suggested that black holes of active (Seyfert 1) galaxies follow a similar relation, 
but apparently with a significantly lower value of AfBH/-A^buigo (Wandel 1999). New data show that this 
difference was mainly due to over-estimating the black hole mass in quiescent galaxies and over-estimating 
the bulge magnitude of Seyfert galaxies. Using new and updated data we show that AGNs (Seyfert 
galaxies and quasars) follow the same BH-bulge relation as ordinary (inactive) galaxies. We derive the 
BH/bulge relation for a sample of 55 AGNs and 35 quiescent galaxies, finding that broad line AGNs 
have an average black hole/bulge mass fraction of ~ 0.0015 with a strong correlation (Mbh oc L^^^^ ^^)- 
This BH-bulge relation is consistent with the BH-bulge relation of quiescent galaxies and much tighter 
than previous results. Narrow line AGNs appear to have a lower ratio, MBH/Afbuigc ~ 10~^ — 10^"^. We 
find this to be a more general feature, the BH/bulge ratio in AGNs being inversely correlated with the 
emission-line width, implying a strong linear relation between the size of the broad emission line region 
and the luminosity of the bulge. Finally, combining AGNs with observed and estimated stellar velocity 
dispersion, we find a significant correlation (A/bh ct"^'^^^), consistent with that of quiescent galaxies. 

Subject headings: black hole physics, galaxies: bulges, galaxies: nuclei, galaxies:active, galaxies: Seyfert, 
quasars:general 



1. INTRODUCTION 

Compact dark masses, probably massive black 
holes (MBHs), have been detected in the cores of many 
nearby galaxies (Kormcndy and Richstone 1995, Kor- 
mendy & Gebhardt (2001)). Magorrian et al. (1998) 
presented a sample of 32 nearby galaxies with MBHs and 
suggested that the black hole (BH) mass is proportional 
to the luminosity of the host galaxy or host bulge (or 
equivalently, to their mass) hereafter referred to as the 
BH/bulge relation, with the BH mass being about 0.006 
of the mass of the spheroidal bulge (although with a sig- 
nificant scatter). 

In addition to the MBHs detected by techniques of stel- 
lar and gas kinematics, the masses of about three dozen 
MBHs in AGNs have been estimated by reverberation 
mapping of the broad emission-line region. High qual- 
ity reverberation data and virial BH mass estimates are 
presently available for 20 Seyfert 1 nuclei (Wandel, Peter- 
son and Malkan 1999, hereafter WPM) and 17 PC quasars 
(Kaspi et al. 2000). Laor (1998) suggested that the masses 
of BHs in quasars (estimated empirically) and their host 
bulges follow the BH/bulge relation of Magorrian et al. 
(1998). Using the WPM Seyfert sample with bulge esti- 
mates from Whittle (1992) Wandel (1999) found that the 
BH/bulge ratio is significantly lower (on average, by a fac- 
tor of ^20) than the ratio found by Magorrian et al. (1998) 
for MBHs in quiescent galaxies. 

Recent research using higher quality HST data and a 
more careful treatment of the modeling uncertainties give 
lower values for black hole masses in nearby galaxies (Mer- 
ritt and Ferrarese 2001a), with an average black hole to 



bulge mass ratio of ^0.001 and a nearly linear BH-bulge 
relation Mbh oc ibuige (Kormendy & Gebhardt 2001). 

The updated lower dynamical BH mass estimates for 
quiescent galaxies reduced the discrepancy between inac- 
tive galaxies and reverberation mapped Seyferts, but Con- 
sequently, the quasars in Laor's (1998) sample appear to 
have larger BH-bulge ratios than inactive galaxies and of 
less luminous AGNs. To explain this Laor (2001) sug- 
gested a nonlinear BH-bulge mass relation, Mbh ~ -^bu^ge- 

Ferrarese and Merritt (2000) and Gebhardt et al. 
(2000a) have found that the mass of MBHs in quiescent 
galaxies is better correlated with the stellar velocity dis- 
persion in the bulge than with the bulge luminosity. Al- 
though to date there are only a few (less than ten) Seyfert 
galaxies with reliable measurements of stellar velocity and 
with reverberation BH mass estimates (Ferrarese et al. 
2001), they seem to be consistent with the BH- velocity 
dispersion relation of inactive galaxies. 

Recently McLure & Dunlop (2001, hereafter MDOl) 
have studied the BH-bulge relation of 30 quasars, com- 
bined with the Seyfert galaxies in the WPM sample for 
which they estimate bulge magnitudes from HST imaging 
and bulge-disk decomposition. Assuming a flattened BLR 
geometry and a biased distribution of inclinations MDOl 
suggest that the actual BH mass of AGNs is 3 times larger 
than inferred from the reverberation virial method assum- 
ing isotropic geometry. While at present a factor of 2-3 in 
the BH reverberation mass cannot be ruled out, system- 
atically increasing the Seyfert BH masses by a factor of 3 
seems to violate the agreement between Seyferts and qui- 
escent galaxies in the BH-velocity dispersion relation (fig. 
2 of Ferrarese et al. 2001). 



2 



Wandel 



In this work we reconcile the BH-bulge mass relation of 
AGN with that of quiescent galaxies, analyzing the fac- 
tors which led to the earlier apparent discrepancies be- 
tween Seyferts, quasars and inactive galaxies, and demon- 
strate that all three groups do follow the same BH-bulge 
mass relation without a spatial geometry or bias for AGNs. 
In section 2 we present the sample, including 23 Seyfert 
galaxies, 32 quasars and 35 inactive galaxies with reliable 
Mbh and bulge data. In section 3 we show that the BH- 
bulge mass relation of the combined Seyfert -|-quasar sam- 
ple is consistent with that of inactive galaxies, both in 
slope and in normalization, discuss the uncertainties and 
point out the factors that produced the apparent discrep- 
ancy in the Mbh — -Mbuigo relation between Seyferts and 
inactive galaxies in previous works. We also find that Nar- 
row Line Seyfert Is and narrow line quasars have signifi- 
cantly lower MeH/Mbuigo ratios than broad line AGNs and 
inactive galaxies, which is the result of a systematic trend 
(section 4) - a strong correlation between the BH mass 
to bulge luminosity (or mass) ratio and the emission line 
width of the active nucleus. In section 5 we analyse the 
implications of the BH mass - stellar velocity dispersion 
relation, while section 6 discusses possible reasons for the 
lower BH/bulge ratio of narrow line AGNs. 



Seyferts with the broadest hues (F9, IC 4329A, NGC 3227 
and NGC 4151) have correction factors of order 10, the 
narrow line objects (MRK 335 and the marginally narrow 
line 3C120) have a negligible correction, « 1. The best fit 
to the bulge luminosity correction for the 15 Seyferts in 
MDOl is 

/ FWHM(H/3) \ 
V 103km s"i / 



l0g(Lsdv/^iing) = 1.3 log 



0.17 (1) 



The correlation coefficient is R=0.77 with a < 0.1% prob- 
ability of chance correlation. This empirical correlation is 
so tight that we use it to predict the actual bulge magni- 
tude from the SdV relation for Seyferts which do not have 
directly measured bulge luminosity. With the empirical 
formula we obtain bulge magnitudes for the five Seyferts 
in our sample which do not have measured bulge lumi- 
nosity (AKN 120, MRK 110, MRK 279, NGC 3516, NGC 
4593). In the above analyses we have used for 3C390.3 
a FWHM H/3=3500km s'^ (Crenshaw 1986) which is sub- 
stantially lower than the value of lO^km s~^ reported in 
WPM; this very broad feature may be temporary, due to 
the line profile variability (Peterson, private communica- 
tion). 



2. DATA 

We compile a sample of 55 AGNs with BH mass and 
bulge magnitude. The AGN data are listed in tables 1 
(Seyferts) and 2 (quasars). We use Ho — 80km s~^Mpc~^ 
and go = 0.5. We convert magnitudes from other works to 
the V band, using the relations V — R = 0.8 and B — V = 
0.8. We express the bulge magnitude in units of Lq us- 
ing the standard relation ibuigo/-^© — 0.4(— My -I- 4.83). 
To find the bulge mass we use the relation M/Mq — 
1.18(L/Lo) - 1.11 (Magorrian et al. 1998). 

2.1. Seyferts 

We use the sample of 17 Seyferts in Wandel, Peter- 
son and Malkan (1999), supplemented by 6 more Seyferts 
with reliable BH mass or bulge magnitude data: MRK 
279 (Santos Lleo et al. 2001), NGC 3516 and NGC 4593 
(Ho 1999, Gebhardt et al. 2000b) and MRK 841, NGC 
4253 and NGC 6814 (bulge magnitudes from Virani et 
al. , BH mass from L-R relation (cq. 2 below, cf. Laor 
2001). For the Seyfert galaxies in the WPM sample we 
use the bulge magnitude given by MDOl, estimated from 
HST imaging and two dimensional bulge-disk decomposi- 
tion. This method is more accurate than the empirical 
formula for the bulge/total ratio depending on the Hubble 
type (Simien & deVaucouleurs 1986; Whittle et al. 1992) 
used in Wandel (1999). The bulge luminosity obtained by 
bulge/disk decomposition of galaxy images is systemati- 
cally lower than the empirical estimate: for the 15 Seyferts 
common to Wandel (1999) and MDOl the latter have an 
average bulge magnitude lower by a factor of 3. However, 
the difference in the bulge magnitude estimated by the two 
methods is not uniform: for individual objects it varies be- 
tween 1 and 10. We find that the correction factor (that 
is, the ratio of the bulge luminosity derived by MDOl from 
HST imaging and that estimated from the Simien & de- 
Vaucouleurs (SdV) relation) is correlated with the FWHM 
of the broad emission lines (fig. 1). For example, while the 
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Fig. 1. — The bulge luminosity correction: the difference be- 
tween the bulge/disc decomposition method (MDOl) and the em- 
pirical Simien-de Vaucouleurs method (bulge/total dependence on 
Hubble type) vs. the H/3 line width for Seyfert 1 galaxies. 

2.2. Quasars 

We adopt the quasar sample of MDOl, supplemented 
by PG1425-H267 (Kirhakos et al. 1999) and PG1704-H60 
(BH mass from Kaspi et al. 2000, host magnitude from 
Hamilton, Casetano and Turnshek 2000). 

The virial formula for the BH mass in AGNs requires 
the size and gas velocity of the broad emission line region. 
Most of the quasars do not have reverberation data. For 
those we estimate the size of the BLR using the empirical 
correlation 

RsLR = 33(i/L^(5100)/10'*Vg s"i)" '^light - days (2) 

between the BLR size and the continuum luminosity 
(Kaspi et al. 2000). We test the goodness of the empiri- 
cal R-L relation for the quasars with reverberation data. 
There are four quasars in the MDOl sample with reverber- 
ation data (Kaspi et al. 2000) - PG0052, PG0953, 3C273 
and PG1307. Comparing the BH masses the agreement 
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Name 


FWHM 


logLv.T 


logLv.img 


log(Mbui) 


log(MBH) 


1 / ^BH \ 


cr(km s ) 


<51ogL|,„i 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


3C 120 


2.21 


9.97" 


10.03 


10.72 


7.49 


-3.23 


162 


0.52 


3C 390.3 


10.5 


10.70 


10.18 


10.90 


8.59 


-2.31 




0.52 


Akn 120 


5.85 


10.28 


9.45 ° 


10.04 


8.29 


-2.65 






F9 


5.90 


10.75 


9.81 


10.46 


7.94 


-2.52 




0.94 


IC 4329A 


5.96 


9.82 


8.86 


9.34 


< 6.86 


< -2.47 




0.97 


Mrk 79 


6.28 


9.93 


9.54 


10.15 


8.02 


-2.13 


125 


0.38 


MrkllO * 


1.67 


10.16° 


9.70° 


10.74 


6.91 


-3.83 


90 




Mrk 279 


3.41 


10.24 


9.72° 


10.36 


7.41 


-2.95 






Mrk 335* 


1.26 


9.86" 


9.89 


10.56 


6.58 


-3.98 




-0.03 


ivirK ouy 




lu.oo 


1 n nfi 

lU.UD 


lu. / D 


7 QS 






u.ou 




2.17 


10.36 


10.29 


11.03 


7.15 


-3.88 


169 


0.07 


Mrk 817 


4.01 


10.32 


9.69 


10.33 


7.56 


-2.77 


140 


0.63 


iVl 1 K. O -L 


5 41 




9 69 


10 91 


8 49^ 


-2 42 






i> v_l V_' 1 


5 53 


10 04 


8 96 


9.47 


7 69 


-1 78 


128 


1.07 


NGC 3516 


2.70 


9.61 


10.00° 


10.69 


7.36 


-3.33 


124 




NGC 3783 


4.10 


9.88 


9.15 


9.69 


7.04 


-2.65 




0.73 


NGC 4051* 


1.23 


9.73 


9.46 


10.05 


6.20 


-3.85 


80 


0.27 


NGC 4151 


5.23 


9.84 


8.83 


9.31 


7.08 


-2.23 


90 


1.01 


NGC 4253 


4.10 




9.63 


10.25 


7.00'' 


-3.25 






NGC 4593 


3.72 


9.59 


10.20° 


10.21 


o.gi*" 


-3.30 


124 




NGC 5548 


5.50 


10.21 


9.68 


10.31 


7.83 


-2.48 


180 


0.53 


NGC 6814 


5.50 


9.56 


9.25 


9.81 


7.08*' 


-2.73 




0.31 


NGC 7469 


3.20 


10.21 


10.04 


10.74 


6.88 


-3.86 




0.17 



Unknown Hubble type, bulge correction estimated assuming Sa 

Bulge magnitude from Virani (1999), BH mass from R-L scaling (cq. 2) 
No image, bulge luminosity estimated from FWHM-5Lbuigo relation (eq. 1) 

Note. — Narrow line Seyferts are denoted by an asterisk. Column (2) - FWHM of (H/3), rms profile, in units of lO^km s~^. (3) - absolute 
bulge V luminosity derived from Hubble Type empirical formula from Whittle et al. (1992), (4) - absolute bulge V luminosity derived by 
imaging bulge/disk decomposition from MDOl, (5) — log of the galactic bulge mass (M(,„;) in Mq (MDOl), (6) - black hole mass (from WPM), 
(7) - BH to bulge mass ratio, (8) - stellar velocity dispersion in km s~^from Nelson & Whittle (l995) and Ferrarese et al. (2001), (9) - log 
bulge correction (Ly^x/ Limg)- 



Table 2 
Quasar data 



Name 


FWHM 


logLv 


log(Mbui) 


log(MBH) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


0052-1-25 


4.37 


10.41 


11.18 


8.29 


-2.89 


0054-1-14 


9.66 


10.65 


11.46 


8.91 


-2.55 


0137-1-012 


7.61 


10.81 


11.65 


8.57 


-3.08 


0157-1-001* 


2.14 


10.93 


11.79 


7.70 


-4.09 


0204-1-292 


6.80 


10.53 


11.32 


7.15 


-4.17 


0205-1-02 


2.90 


9.65 


10.28 


7.86 


-2.42 


0244-1-194 


3.70 


10.33 


11.08 


8.03 


-3.05 


0736-1-017 


2.96 


10.65 


11.46 


7.99 


-3.47 


0923-1-20 


7.30 


10.53 


11.32 


8.94 


-2.38 


0953-1-41 


2.96 


10.33 


11.08 


8.39 


-2.69 


1004-1-13 


6.34 


10.85 


11.70 


9.09 


-2.61 


1012-1-008 


2.64 


10.77 


11.60 


7.79 


-3.81 


1020-103 


7.95 


10.57 


11.36 


8.35 


-3.01 


1029-14 


7.50 


10.41 


11.18 


9.07 


-2.11 


1116-1-215 


2.92 


10.69 


11.51 


8.22 


-3.29 


1202-1-28 


5.01 


10.33 


11.08 


8.29 


-2.79 


1217-1-023 


3.83 


10.73 


11.55 


8.40 


-3.15 


1226-1-02 


3.52 


10.97 


11.84 


8.61 


-3.23 


1302-10 


3.40 


10.61 


11.41 


8.30 


-3.11 


1307-1-08 


5.32 


10.37 


11.13 


8.05 


-3.08 


1309-1-35 


2.94 


10.41 


11.18 


7.99 


-3.19 


1402-1-26* 


1.91 


9.93 


10.61 


7.28 


-3.33 


1425-1-267" 


9.40 


10.57 


11.36 


9.32 


-2.04 


1444-1-40 


2.48 


10.37 


11.13 


8.06 


-3.07 


1545-1-21 


7.03 


10.53 


11.32 


8.93 


-2.39 


1635-1-119 


5.10 


10.49 


11.27 


8.10 


-3.17 


1704-1-60*'' 


0.40 


10.37 


11.13 


6.87 


-4.86 


2135-14 


5.50 


10.61 


11.41 


8.94 


-2.47 


2141-1-175 


4.45 


10.73 


11.55 


8.74 


-2.81 


2247-1-140* 


2.22 


10.73 


11.55 


7.59 


-3.96 


2349-01 


5.50 


10.97 


11.84 


8.78 


-3.06 


2355-082 


7.51 


10.65 


11.46 


8.39 


-3.07 



Kirhakos et al. 1999 

FWHM(rms) and M(rms) from Kaspi et al. (2000) 



Note. — Narrow-line quasars are denoted by an asterisk. Column (2) - FWHM of (H/3) in units of lO^km s ^. (3) - absolute bulge V 
luminosity (in units of L0,from MDOl), (4) - log of the galactic bulge mass, (5) - h\3x:k hole mass, (6) - BH to bulge mass ratio 
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between the empirical R-L relation and the reverbera- 
tion result is within a factor of 1.5. 

There are three quasars for which our BH mass (calcu- 
lated with the empirical formula) differs significantly from 
the values given by MDOl (after correcting for their differ- 
ent kinematic factor, see 2.4 below): PG1307 (the mass in 
MDOl is by a factor 4 too low, apparently because a mis- 
taken FWHM), similarly for 0204-^292 (3C059) MDOl get 
a BH mass too low by a factor 40 (they have FWHM=1040 
km/s, while we use a FWHM(H/3) of 6800 km/s, cf. Er- 
ackleous & Halpern 1994) and PG 1012-1-008 (their BH 
mass is by a factor 10 too low, probably a typo). 

2.3. Quiescent Galaxies 

For comparison with inactive or weakly active galaxies 
we adopt the sample of Kormendy & Gebhardt (2001), 
which lists 37 galaxies with BH mass determinations from 
stellar or gas dynamics and bulge magnitudes. We use the 
relation B-V=0.8 in order to translate the B bulge mag- 
nitude to our standard V luminosity. We exclude the two 
galaxies NGC4486B and NGC5845 which have a large un- 
certainty (factor 10) in the BH mass, significantly larger 
than the other objects in the sample. In the Mbh — c 
analyses we also exclude NGC 4945 which does not have 
a a observation. 

2.4. Black Hole Mass 

The BH mass of AGNs can be estimated assuming the 
emission-line gas is bound (see section 6.5). The virial re- 
lation the gives M = k'^G~^RbirVFWHM^ where kvpwHM 
is the velocity dispersion in the gas (deduced from the ob- 
served width vpwHM of the emission lines) and the kine- 
matic factor fc, which depends on the kinematics, geometry 
and emissivity of the BLR. Assuming the velocity disper- 
sion in the line emitting gas is isotropic we have k = 
which gives 

A/bh = 1.46Me X 10^ ( ( V (3) 

° Vlt - days; V 103km s^V 

Non isotropic geometry and velocity dispersion (combined 
with a non-uniform distribution of inclinations) can give 
different values of the coefficient k (see section 6.4). MDOl 
assume a flattened BLR geometry and a distribution of in- 
clination biased to face-on, which yield fc = 3/2 and hence 
the virial BH mass is 3 times larger than in the isotropic 
geometry. 

3. THE BH-BULGE RELATION 

3.1. The Luminosity Relation 

In fig. 2 we plot the BH mass vs. bulge luminos- 
ity. Notably the scatter in the Seyfert and quasar groups 
is quite large. We note however that all AGNs with 
low MBn/^buigc values are narrow emission- line objects 
(NLSls or narrow line quasars). In other words, it seems 
that narrow line AGNs tend to have small BHs compared 
to their host bulges. We draw the limit between narrow 
and broad line Seyferts at 2200km s^^, which is close to 
the traditional NLSl definition (2000km s^^) but includes 
a couple of narrow line objects in our sample which are 
very close to that definition (MRK 590, PG0157-H00 and 
PG2247-I-14). The precise value is not important, as the 



NLSl definition is merely a historic convention, and our 
results are not very sensitive to the exact border value. 

In section 4 we show that there is a strong correlation 
between line width and BH/bulge ratio, so that the tran- 
sition from broad to narrow line objects is gradual, as well 
as the transition from high to low BH/bulge ratios. This 
correlation further supports the separation of the narrow- 
line objects: as the BH mass estimate of AGNs strongly 
depends on the line width, removing the narrow-line ob- 
jects may seem to automatically remove the outlying data 
points with low BH/bulge ratio. However, we find that 
this is not the the BH/bulge mass ratio depends 

also on other observables: the BLR radius and the bulge 
magnitude. In order to test the significance of separating 
the narrow line objects, we compare the two correlations: 
BH mass vs. FWHM and BH/bulge ratio vs. FWHM. 
The correlation coefficients are 0.62 and 0.80, respectively, 
which shows that the BH/bulge - FWHM relation is more 
significant than the mere BH mass dependence on the 
FWHM. 

In section 6 we discuss a few physical and observational 
models which could explain such a correlation and moti- 
vate the separation of the narrow line objects. We group 
those objects (Narrow Line Seyfert Is and narrow line 
quasars) as a separate class which we denote NLAN (Nar- 
row Line Active Nuclei). 
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Fig. 2. — The BH mass vs. bulge luminosity. Green diamonds de- 
note BHs detected by stellar and gas dynamics in quiescent galaxies 
(Kormendy & Gebhardt 2001), blue circles - quasars (MDOl), red 
triangles - Seyfert 1 galaxies (BH mass from WPM, bulge luminosity 
from MDOl) and pink squares - Narrow Line Active Nuclei (NLSls 
from WPM and NLQs from MDOl). The least-squares linear fit for 
quiescent galaxies is indicated by a dashed green line and NLANs 
by a pink, dot dashed line. The arrows indicate the effects of the 
improvements in the data that cancelled the apparent former dis- 
crepancy in the Mgjj/L^^^i^^ ratio between Seyferts and quiescent 
galaxies and quasars (see text). 

The best fits to the data are (i and M are in solar units) 

1. Quiescent galaxies (35): 

log A/ = (0.96 ± 0.13) logi- (1.7 ± 1.3), 

2. Broad line AGNs (19 Seyferts and 28 quasars) 
log A/ = (0.90 ± 0.11) logL- (1.1 ± 1.1), 
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3. Narrow line AGNs (4 Seyferts and 4 quasars): 
logM = (0.84 ± 0.21) logL- (1.5 ±2.2) 

4. Quiescent galaxies and broad line AGNs (82) 
log A/ = (0.93 ±0.13) log L- (2.9 ±2.6) . 

The fits for the BH/bulge mass relation, along with the 
corresponding correlation measures are detailed in table 3 
below. The (broad line) AGNs have an almost linear fit 
very similar to that of inactive galaxies and to the com- 
bined sample of broad-line AGNs and quiescent galaxies. 
The NLAN group has a similar slope but with an M/L 
ratio approximately a factor of 10 lower. 

3.2. The Mass Relation 

Transforming bulge luminosity to mass one finds a sim- 
ilar correlation between the BH mass and the bulge mass. 
The correlation is analogous to the BH mass vs. bulge 
luminosity relation (fig. 2), with a somewhat less steep 
slope (since Mbuigc ~ L^'^). Fig. 3 shows the BH mass 
vs. the bulge mass. In we list the best fits to the A/bh 
vs. Mbuigo relation for broad line AGNs (Seyfert Is and 
quasars), narrow line AGNs (NLSls and NL QSOs) and 
quiescent galaxies, as well as for the combined sample. For 
each subset we give the three measures of the fit: the corre- 
lation coefficient, the standard deviation of the (logarith- 
mic) scatter, given by dM = [Y.{5\og{Mbhf / {N - 1)]!/^ 
and the reduced of the fit. Since the BH mass esti- 
mates have systematic modeling uncertainties (section 4) , 
we approximate the uncertainty in the calculation by 
a uniform error of 0.5 dex (factor 3) in the BH mass es- 
timate, which is derived by combining the typical obser- 
vational error with the systematic uncertainty (see section 
3.3 below). 

We also list the mean BH/bulge mass ratios and their 
standard deviations. In order to compare the BH-bulge 
correlation with that of the Afen — cr relation we have 
added a row showing the A/bh — cr correlation for qui- 
escent galaxies. Since both correlations have the same 
y- variable (A^bh) they can be compared using the same 
statitstic. For consistency we treat the errors in the BH 
mass for AGNs and quiescent galaxies in the same manner, 
although in some cases, such as the Galaxy and NGC 4258, 
the error in the BH mass is much smaller. Note that the 
effective error (factor 3) assumed for the BH mass is larger 
than the measurement uncertainty for most BHs in quies- 
cent galaxies. Using merely the measurement uncertainty 
in A/bh, e.g. the error bars quoted in the Kormendy and 
Gebhardt (2001) sample gives =3.5 for the A^bh — cr 
relation and x^ — 10.5 for AfBH-bulge relation. 

We note several important results: 

1. Broad- line (BL) AGNs have a very similar 
distribution (for both, slope and mean MH/bulge 
ratio) to that of the quiescent galaxies. 

2. In spite of their small number in out sample, 
removing the narrow line AGNs has a significant 
effect on the combined AGN sample: the total 
AGN sample has a significantly lower correlation 
coefficient, and a lower mean BH/bulge ratio than 
BL AGNs and quiescent galaxies. 



3. NLANs have an average a BH/bulge ratio lower by 
a factor of 10 than the broad line AGN average, 
and the normalization of the best fit is similarly 
lower. 

4. Quasars and Seyfert galaxies as separate groups 
have shallower slopes and lower correlation 
coefficients, due to the lower dynamical range of 
each of these groups separately. 

5. The BH-bulge correlation of broad line AGNs 
is quite tight, not much less than the A/bh — cr 
relation for quiescent galaxies. 

The last item is demonstrated by comparing the three 
correlation statistics in table 3 (the correlation coefficient, 
the logarithmic scatter and x^) for the BH-bulge and the 
BH-CT relations. For the 47 BL AGNs we find 0.76, 0.43 
and 0.78, respectively, nearly as tight as the A/bh — f cor- 
relation of the 34 quiescent galaxies in the Kormendy and 
Gebhardt (2001) sample (0.86, 0.38 and 0.60). We note 
that excluding the Seyfert galaxy NGC 7459 (which has 
an exceptionally low BH mass relative to its luminosity) 
the BH-bulgc correlation of the broad-line AGNs becomes 
even tighter (dM==0.39). 
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Fig. 3. — Mass estimates of MBHs plotted against the mass of the 
host bulge (or host galaxy for elliptical galaxies). Green Diamonds 
denote BHs in inactive galaxies, blue circles - quasars, red triangles 
- Seyfert 1 galaxies and pink squares - Narrow Line Active Nuclei 
(NLSls and NLQs). The least-squares linear fits for each class are 
indicated by a solid green line (quiescent galaxies), dashed red line 
(broad line AGNs) and a pink, long dashed line (NLANs). 

3.3. Uncertainties 

The formal uncertainty in the slope of the linear fit to 
the (broad line) AGNs is 0.09. However, for a realistic 
uncertainty it is necessary to take into account the errors 
in the BH mass estimate and bulge luminosity. The er- 
rors in the BH mass are of two kinds: systematic errors 
in the reverberation mapping and virial mass calculation 
and measurement errors. The latter are random and typ- 
ically of the order of 2-3 (see WPM for the observational 
uncertainty in the reverberation masses and for the error 
of the L — R scaling method (used in our sample for most 
of the quasars) compared to the reverberation technique ). 
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Table 3 
BH-BULGE Mass Correlation's 



oet 


N 


a 


b 


MBH/-Mblg 


R 


6A4 


,^2 

X 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


All AGNs 


55 


0.74±0.11 


-0.09±1.22 


-2.93±0.57 


0.67 


0.54 


1.20 


Broad Line AGNs 


47 


0.73±0.09 


-0.15±1.02 


-2.80±0.48 


0.76 


0.43 


0.78 


Quiescent Galaxies 


35 


0.82±0.11 


-0.82±1.15 


-2.77±0.50 


0.80 


0.48 


0.96 


BL AGNs and Q. Galaxies 


82 


0.79±0.07 


-0.53±0.74 


-2.79±0.47 


0.79 


0.45 


0.82 


BL quasars 


28 


0.51 ±0.24 


2.70±2.74 


-2.86±0.41 


0.38 


0.38 


0.61 


BL Seyfert Is 


19 


0.52±0.23 


2.27±2.37 


-2.68±0.55 


0.48 


0.49 


1.01 


NLANs 


8 


0.78±0.17 


-1.53±1.9 


-3.9±0.27 


0.88 


0.24 


0.27 


Quiescent Galaxies-(7 


34 


3.65±0.37 


-0.29±0.85 


5.71±0.60 


0.86 


0.38 


0.60 



Note. — Columns: (1) N - number of objects, (2-3) coefficients of the linear fit log(MBH) = <J log{A^bulge) + b ^■iid standard deviations, (4) 
mean log(BH/bulge) mass ratio and standard deviation, (5) R - correlation coefficient, (6) Standard deviation scatter of the fit i5(log(MBH)) 
(7) - reduced assuming the error in log(MBH) is 0.5 dex. 



The systematic errors are not random, but depend on the 
geometry and dynamics of the broad emission hne region, 
and in extreme cases could amount to a factor of 3 in each 
direction (e.g. Krolik 2001; see section 6 below). However, 
since this dependence is largely unknown and can act in 
both directions, depending on the conditions in each AGN, 
we treat the systematic errors as random ones as well. 

The uncertainty of a linear Ht y = ax + b for N data pairs 
{x, y) with typical errors 6x and 6y may be estimated as 
6a = [^a^^^i + {a6x/y/N)^ + {6y/VNfYI^. Combining 
the formal error in the fit with the random measurement 
errors in the BH mass (assumed to be a factor of 3) and 
bulge luminosity (factor 2), and the systematic errors in 
the BH mass estimate (assumed to be random with a fac- 
tor of 3) for the broad-line AGN sample (N=47), gives 
5a = 0.16, so that 

A^BH « (4a) 

and 

MbH <X M0,^^±0.15 (4^) 

with < MBH/Mbuigo >= 0.0015. 

Studying a sample of 16 AGNs Laor (2001) finds a 
steeper fit, Mbh oc -^bu^g^'^'^'' (formal error only), which 
is ^ 5 — 6(T higher than what we find in this work (A'/bh 
M°J3±o.o9^ This difference is explained as follows, (i) 
Laor's fit includs quiescent galaxies, some of which have a 
low bulge luminosity and a low MBn/A^buigc value (see sec- 
tion 6.2). For the AGNs separately Laor finds a less steep 
slope of 1.36±0.15. Furthermore, a single point signifi- 
cantly influences the slope - the NLSl NGC 4051 - which 
extends the dynamical range by almost an order of mag- 
nitude. If it were excluded (being a narrow line object), 
the slope of the fit would decrease considerably (and the 
standard deviation in the slope would increase). Exclud- 
ing NGC 4051 and adding the measurement error would 



flatten the slope and increase the uncertainty in Laor's fit, 
which would ecome consistent with linearity. 

(ii) The BH masses used by Laor (1998; 2001) arc sys- 
tematically larger than those used in this work (and in 
MDOl). The BH masses in Laor's sample arc arc on aver- 
age larger by a factor of 1.5 than those calculated for the 
same quasars in our sample. This discrepancy originates 
in part from the different luminosity measures used: Laor 
used the broad band (0.1-1 micron) luminosity (Neuge- 
bauer et. 1979) while we (and MDOl) use the monochro- 
matic luminosity j/JV(5100). The latter seems more ade- 
quate here, as we (as well as MDOl) estimate the BH mass 
of the quasars using the empirical R — L relation of Kaspi 
et al. (2000) which has been derived with vF^iblOQ). 

We conclude that quasars, Seyfert galaxies and inactive 
galaxies have the same BH-bulge relation and BH/bulge 
ratio, while NLANs seem to have BH/bulge ratios lower 
by factors of ~l-30. 

3.4. Resolving the discrepancy: Do Seyferts have a lower 
BH/bulge ratio than Quiescent Galaxies? 

Previous works have found a signiflcant difference in the 
black hole mass- bulge luminosity (or mass) relationship 
of Seyferts and ordinary galaxies, the former showing sys- 
tematically lower MBn/^buige values (Wandel 1999, Ho 
1999, Gebhardt 2000b). The discrepancy between the 
MeH/A/buige ratio of quiescent galaxies and Seyferts (fac- 
tor of 10-20 on average) was a result of three factors, shown 
schematically by the arrows in fig. 2: (i) the Magorrian 
BH masses were overestimated by an average factor of ~ 5 
(Men/Afbuige =0.6% compared with the current estimate 
of 0.1-0.13% (Ho 1999; Merritt & Ferrarese 2001a). The 
average decrease in the quiescent galaxy BH mass is shown 
in fig. 2 by the arrow denoted G. (ii) the Seyfert bulge lu- 
minosity calculated using the Simien de Vaucouleurs for- 
mula is too large (by a factor of 2 — 3); the bulge lumi- 
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nosity estimated by direct imaging and bulge/disk decom- 
position is lower by an average factor of ~2 (O.lSdex in 
our total AGN sample, as seen in table 3, but more when 
applied to smaller sub-samples, e.g. to the Seyferts) repre- 
sented by the horizontal arrows, (iii) Narrow Line Seyfert 
1 galaxies, which do seem to have BH/bulge mass ratios 
lower than broad line active galaxies, lowered the aver- 
age. Excluding the NLSls raises the average MBH/A^buigo 
value of AGNs by a factor of '--^ 2, denoted in fig. 2 by 
the arrow pointing upwards. Items (ii-iii) combined in- 
creased the BH/bulge ratio of Seyferts by an average fac- 
tor of ^--^3.5. (dashed arrow). Similarly, the disagreement 
between the BH/bulge ratios of the PG quasars and rever- 
beration Seyferts (Wandel 1999) disappeared due to the 
revision of the Seyferts' BH/bulge ratio and the lower BH 
mass estimates of the quasars (the arrow denoted by Q in 
fig 2). 

Given the present good agreement of AGNs and qui- 
escent galaxies, one may ask whether there is room for 
systematically higher BH mass estimates in AGNs, e.g. 
due to inclination effects. In the geometry suggested by 
MDOl, the virial BH mass of AGN would be larger on av- 
erage factor of 3, compared with the mass calculated under 
the assumption of isotropic velocity distribution. We can 
check what would be the implication of such an increase 
on the BH-bulge and BH-velocity dispersion relations oi 
AGNs compared with those of quiescent galaxies (where 
the BH mass estimates are not affected). In order to test 
this hypothesis we calculate the of the AGN sample 
(including the NLANs) with the BH mass increased by a 
uniform factor 3, with respect to the best fit to the quies- 
cent galaxies presented in table 3. Assuming the typical 
measurement + modeling error in the BH mass is 0.5 dex, 
this gives a reduced of 2.2 compared with 1.20 for BH 
masses calculated with an " isotropic BLR" assumption for 
all AGNs in our sample. A similar exercise can be done 
for the a — Mbh relation, which has less AGNs, but also 
lower spread and a smaller error in the BH mass. Increas- 
ing each BH mass by a factor of 3, we find for the 11 
Seyfert galaxies with velocity dispersion data (table 1) 
of 2.58 compared with 1.08 for the "isotropic BLR" BH 
masses. We conclude that the BH mass calculated under 
the "isotropic BLR" assumption for the BLR geometry is 
more consistent with the BH-bulge and BH-velocity dis- 
persion relations for quiescent galaxies. 

4. THE BH-BULGE VS. EMISSION-LINE WIDTH 
CORRELATION 

We find that narrow-line AGNs - NLSls and narrow line 
quasars - tend to have low MBH/-Mbuigc values. This gives 
a physical motivation to the exclusion of the narrow line 
AGNs as a separate group in the BH-bulge relation. 

Mathur et al. (2001) make a similar suggestion based on 
a sample of 15 NLSls, but with less reliable BH masses. 
Only three of the objects in their sample have black 
hole masses measured by reverberation , 3 more by the 
scaling-virial method, and the 9 remaining BH masses are 
calculated using accretion disk modeling, which is less re- 
liable. They find that their accretion disk model gives BH 
masses too large (by a factor of 3-7) compared with the 
reverberation masses (where available). In order to com- 
pensate for this discrepancy they decrease the BH masses 



given by the accretion-disk model by a calibration correc- 
tion factor of 5. For the bulge magnitude Mathur et al. 
use Whittle (1992) with the Simien-deVaucouleurs rela- 
tion (which is likely to overestimate the bulge magnitude 
in AGN, as discussed in above). 

Considering all 55 AGNs in our sample with avail- 
able -A/bh and bulge magnitude we find that the smaller 
BH/bulge ratio of narrow line objects is in fact the lower 
end of a more general relation: the BH mass /bulge lu- 
minosity ratio is strongly correlated with the emission- 
line width. 
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Fig. 4. — BH-to-bulge luminosity ratio plotted against the H/3 
broad emission-line width for AGNs. Red Triangles represent Seyfert 
1 galaxies and blue circles are quasars. The solid line indicates the 
least-squares linear fit. 

In fig. 4 we plot M^YilU^v\^c vs. FWHM(H/3). The best 
linear least-squares fit is 

V-t^buigo/Q V 10-^km s / 

with a correlation coefficient R — 0.80. BH mass and 
bulge luminosity are in solar units, and we have excluded 
PG1704 because of its highly uncertain reverberation BH 
mass determination (however, if we take the BH mass and 
FWHM values quoted in Kaspi et al. (2000) it agrees well 
with the fit in fig. 4). 

The tight correlation (5) is not surprising - it 
is actually an artifact of the virial relation: since 
A/bh v{FWHM)^ , this correlation merely means that 
RblrJ Lhnigfi is independent on the line width. Combining 
eqs. (3) and (5) we have 

i?BL««27(^^) It -days. 

Looking at the data we find almost precisely the same re- 
sult - a very tight correlation between bulge luminosity 
and BLR size (in preparation), as shown in fig. 5. 

The best fit to all 55 AGNs in our sample combined is 
log(i?/lt - days) = l.O51og(ibuigc/i0) - 8.92 or 

R = 13.5i}o°^ It - days (6) 
with a correlation coefficient of 0.91. We note that while 
Seyferts and NLSls show a significant correlation, quasars 
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alone do not. This is not surprising, as the BLR size of 
most quasars in our sample is determined from the R- 
L scaling relation, which is less reliable than reverbera- 
tion mapping, used for most of the Seyferts. However, the 
quasars, which do have reverberation measured BLR sizes 
(solid circles in fig. 5) are indeed very close to the fit. 
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Fig. 5. — The size of the broad emission-Une region vs. the bulge 
luminosity. Red Triangles are Seyfert 1 galaxies, blue circles are 
quasars and pink squares are narrow line AGNs. Solid symbols indi- 
cate BH masses calculated from reverberation -mapping data (from 
WPM and Kaspi et al. 2000) and open symbols - calculated from 
the L-R empirical relation. The solid line indicates the least-squares 
linear fit. 

This is a non-trivial correlation, not just a derivation, 
as it relates two independent observables: on the one hand 
the bulge luminosity, a global galactic property on a kpc- 
scale, and on the other hand the BLR size (measured by 
reverberation mapping or luminosity scaling), on a few 
light-days scale. 

It may be reflecting the BH-bulge relation, but the cor- 
relation being much stronger than the BH-bulge relation 
supports the case that this new relation is more fundamen- 
tal. A possible connection to basic physical properties may 
be implied by combining the new BLR-bulge relation, the 
empirical relation between the BLR size and the central 
source luminosity (e.g. Kaspi et al. 2000) and the BH- 
bulge relation giving that the Eddington ratio increases 
with BH mass. 

5. THE BLACK HOLE MASS - BULGE VELOCITY 
DISPERSION RELATION FOR AGN 

We have shown that the black hole /bulge ratio in NLSl 
galaxies seems to be significantly smaller than in broad line 
AGNs and in quiescent galaxies. However, Seyfert galaxies 
seem to be consistent with the a — Mbh relation of inac- 
tive galaxies (Gebhardt et al. 2000b; Merritt & Ferrarese 
2001b; Ferrarese et al. 2001). In order to unveil the origin 
of the lower BH/bulge ratio of NLANs (namely, whether 
they really have lower BH masses or rather larger bulges, 
relative to the broad line AGN) we look at their location 
in the a — A/bh plane. The bulge velocity dispersion is 
measured only for 3 NLSls (Mrk 110, Mrk 590 and NGC 
4051) and 8 broad-hne Seyferts (table 1). 

We may use the Faber- Jackson relation to estimate the 
velocity dispersion of other AGNs with measured or es- 



timated bulge luminosity. Also the relation between the 
the narrow emission lines and the velocity dispersion (Nel- 
son 2000) may be used to estimate the velocity dispersion 
in the bulge. Here we use the standard F-J relation, e.g. 
L = iocl, where a2 = cr/lOOkm s~^ and Lq is a luminos- 
ity coefficient determined by a linear fit with a slope of 4 
to the observed Seyfert Galaxies (fig. 6). 

In our limited data set we note that NLSls have a signif- 
icantly larger luminosity coefficient than broad line ones: 
while for broad line Seyfert 1 galaxies Lo = lO^erg s~^, for 
NLSls galaxies Lo = 6 x lO^erg s~^. In addition, galaxies 
with a massive BHs and AGNs in particular, may have 
a different (flatter) F-J relation than galaxies in general 
(Wandel 2001 and work in preparation). 

Using these relations one may compile a larger sam- 
ple and increase the dynamical range by including also 
quasars, to estimate the a — M-qu relation for AGNs. In 
flg. 7 the BH mass is plotted against the stellar veloc- 
ity dispersion. We see that AGNs are consistent with the 
M-QYi — a relation of inactive galaxies (Mbh « cr", with 
a = 3.5-5; Gebhardt et al. (2000a) find a =3.65 and Mer- 
ritt and Ferrarese (2001a) give 4.72). Adding the quasars 
seem to favor a steeper slope, but this result merely re- 
flects the assumed L oc a* relation; if, for example, a flat- 
ter cr — L relation is assumed, the luminous objects would 
have a larger estimated velocity dispersion yielding a flat- 
ter slope for the a — Mbh relation. In flg. 7 we indicate 
by an arrow how the velocity dispersion estimate of lumi- 
nous objects would move if the best flt to the F-J relation 
{L ~ a^-^) for the Seyferts (dotted line in flg 6) were used. 




Fig. 6. — Bulge luminosity plotted against stellar velocity disper- 
sion for Seyfert galaxies. Blue Triangles indicate broad-line Seyfert 
Is, pink squares denote NLSls. The solid lines denote the Faber- 
Jackson relation (L ~ ct*) normalized for the broad and narrow line 
Seyferts, respectively, while the dashed line is the best fit to the 
broad line Seyferts. 



In agreement with previous results (Gebhardt 2000b, 
Merritt & Ferrarese 2001b) the NLANs (sohd and open 
pink squares in flg 7) are more or less consistent with the 
Mbh — cr relation, although they do lie at the lower end of 
the distribution. 

The agreement of the A/bh — f relations of inactive 
galaxies and AGNs (broad and narrow line) suggests that 
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the virial masses are essentially correct also for NLANs. 
The tendency of the NLANs to have low Mbh/c ratios 
suggests that their lower BH-bulge ratios go along with rel- 
atively lower BH masses than broad line AGNs. In order 
to settle this issue, more stellar velocity-dispersion mea- 
surements of AGNs are required, in particular of NLSls. 



10 



9 



7 



6 











MF01 / 












o 

o 






^^GOO 




□ 


o 

/ o 
/ ^6 

^ y o 


D 

O 






□ 


y 

♦ 


o — ^ — 

♦ 














□ 


■ 




♦ BL Sy1 

□ NLAN 




o „Py- 




i ♦ □ 






o BLSy(L) 
o quasars 






□ 















1.8 2.0 2.2 2.4 2.6 



log V* (km/s) 

Fig. 7. — Black hole mass of AGNs plotted against the stel- 
lar velocity dispersion (v* = cr). Blue triangles are broad line 
Seyferts, pink squares denote NLSls and red circles denote quasars. 
Solid symbols dennote Seyferts with measured cr, open ones denote 
Seyferts for which a has been estimated from the Faber-Jackson 
relation (see text). The dashed and solid lines show the u — Mbh 
relation of galaxies (Gebhardt et al. 2000a and Ferrarese & Merritt 
2000, respectively). The arrow shows how the velocity dispersion 
estimate of luminous objects would typically move if the best fit to 
the F-J relation for the Seyferts (fig 6) were used. 

The difference between NLANs and broad line AGNs 
in the BH-bulge relation is much larger than the differ- 
ence in the BH- velocity-dispersion relation. This may be 
related to the bulge mass (and luminosity) being integral, 
extended properties, while the stellar velocity dispersion is 
more closely related to the inner part of the bulge. There 
is more evidence that the black hole mass in AGN is cor- 
related with the velocity dispersion (and hence the virial 
mass) on the scale of the central part of the bulge - the 
narrow line region. This is supported by the correlation 
between the black hole mass and the velocity dispersion 
of the narrow emission-line g measured by the [OHI] 
lines in Seyfert 1 nuclei (Wandel and Mushotzky 1986; 
Nelson 2000). 

6. DO NLANS REALLY HAVE LOWER BH/bULGE RATIOS? 

We have seen that broad line AGNs follow a similar BH- 
bulge relation as inactive galaxies. The question remains 
what can be the cause of the lower BH/bulge ratio in the 
narrow- line AGNs, and whether this effect is real or ap- 
parent. 

Several effects have been mentioned that could cause 
such an apparent lower BH/bulge ratio (e.g. Wandel 
2000): (i) the reverberation - virial method may system- 
atically underestimate the BH mass in NLANs, (ii) the 
measured bulge luminosity in NLANs may be systemati- 
cally too large. Alternatively, the difference between the 
two groups could be real: either (iii) NLANs do have in- 



trinsically larger (or brighter) bulges, or (iv) they have sys- 
tematically smaller black holes . In section 5 we have seen 
that effects (ii) and (iii) are supported by the a — Mbh re- 
lation, as the narrow-line Seyferts with measured velocity 
dispersion (MRK 110, MRK 590 and NGC 4051) seem to 
have larger bulge luminosity than broad line Seyferts with 
comparable velocity dispersions (fig. 6). On the other 
hand, also (i) and (iv) cannot be ruled out, as the rever- 
beration virial BH masses of MRK 590 and NGC 4051 are 
by a factor 2-3 lower than the average Mbh — cr relation 
of broad line Seyferts and inactive galaxies (fig. 7). 

There are a number of observational errors that may 
cause an apparent difference in the black hole - bulge lu- 
minosity relation between NLANs and classic AGNs and 
inactive galaxies, which we discuss below. 

6.1. Errors in estimating bulge luminosity 

As suggested above, the lower Mbh — ^buigo values of 
NLAN could be accounted for if their bulge luminosity are 
systematically over estimated. As we saw, this was indeed 
the case when the Simien-deVaucouleurs method was used 
for Seyfert galaxies. Could a similar problem exist in the 
NLAN set? Two of the 3 NLSls with measured veloc- 
ity dispersion have low Mbh/c values which may indicate 
they have relatively small black holes . However, the loca- 
tion of the NLANs in the Mbh — cr plot (fig 7) with respect 
to the inactive and broad line AGNs indicates that the BH 
mass can account for only a small part of the difference, 
suggesting that NLSls tend to have relatively large bulges. 
It is however difficult to see why the bulge magnitude of 
NLANs be systematically over-estimated by a factor of ~10 
(required to explain the difference in the BH/bulge ratio) 
compared with the broad line AGNs in the sample. 

6.2. Bias introduced by the stellar kinematics method 

Is it possible that NLANs represent a larger population 
of galaxies with low black hole mass to bulge luminosity ra- 
tios, which is under-represented in the presently available 
MBH sample? This could be the case for MBHs detected 
using stellar dynamics, because this method cannot detect 
small black holes , and the detection-limit increases with 
distance. In a resolution-limited method, this would infer 
a lower detection limit which increases with luminosity. 
This hypothesis is supported by the distribution of dynam- 
ically estimated black hole masses in quiescent galaxies. 
There are only four inactive galaxies with BH masses un- 
der IO^Mq: the Milky Way, M32 and NGC 7457 detected 
using stellar dynamics and NGC4945 - with maser dynam- 
ics. The two latter galaxies also have low black hole -to 
bulge luminosity ratios, comparable to the NLAN average 
ratio. In angular-resolution limited methods, the MBH de- 
tection limit is correlated with bulge luminosity: for more 
luminous bulges the detection limit is higher, because the 
stellar velocity dispersion is higher (the Faber-Jackson re- 
lation). In order to detect the dynamic effect of a MBH 
it is necessary to observe closer to the center, while the 
most luminous galaxies tend to be at larger distances, so 
for a given angular resolution the MBH detection limit is 
higher. This may imply that the sample is biased towards 
larger MBHs, as present stellar-dynamical methods are in- 
effective for detecting MBHs below ^ lO^M© (except in 
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the nearest galaxies). Since the BLR method is not sub- 
ject to this constraint, NLANs may represent a low- mass 
BH population. 

6.3. The black hole influence on the stellar velocity 
dispersion 

In section 5 we have seen that the NLANs are not sig- 
nificantly outlying the a — Mbh relation, while they are 
located far below the BH-bulge relation. In principle this 
could be caused by the BH influence on the measured ve- 
locity dispersion. In NLANs the region influenced by the 
BH would be relatively smaller (either because the BH has 
a lower mass or because the bulge is larger, and hence has a 
larger velocity dispersion). This can reduce the difference 
between narrow and broad line AGNs when a rather than 
-t-buige is considered. However, as we show below, the BH 
influence turns out to be negligible for almost all AGNs 
in our sample. The effect of the massive BH on the stel- 
lar velocity dispersion can be estimated by comparing the 
expected velocity fleld (v oc {MBn/rY^^) due to the MBH 
and the velocity dispersion of the host bulge. The black 
hole enhances the velocity dispersion in an observed re- 
gion if its mass is bigger or comparable to the stellar mass 
within the radius corresponding to the projected angular 
size. This can be measured by the dimensionless quantity 

where a is the stellar velocity dispersion in the bulge of 
the host galaxy, 9 is the angular size (e.g. the width of 
the slit) and D is the distance to the galaxy. In other 
words, m is roughly the ratio of the BH mass to the stel- 
lar bulge mass inside the radius corresponding to the an- 
gular size being sampled. The nuclear velocity dispersion 
is typically sampled with a slit width of l-2"(Ferrarese et 
al. 2001 use a slit of 2x4"). Since the measured velocity 
dispersion is weighted by the brightness along the linc-of- 
sight, and since the brightness-density increases steeply at 
the center, the efi'ective value of 6 is probably smaller than 
the slit width. The BH influence parameter can be writ- 
ten as m = 0.9{M/10^MQ)a-'^{D/10Mpc)-'^, where has 
been assumed to be 1". We have calculated m for the 11 
Seyferts with velocity dispersion measurements. The three 
NLSls have m ~ 0.01, lower than the BL Seyferts which 
have typically m ~ 0.03 — 0.1, with the exception of NGC 
3227 which has m 0.24. We conclude that the BH in- 
fluence on the velocity dispersion in the observed AGNs is 
negligible. 

6.4. Inclination and BLR geometry 

If the broad emission-line region (BLR) has a flattened 
geometry and the distribution of inclination angles is bi- 
ased towards face-on inclinations (as one might expect in 
the unified model) the velocity inferred from the observed 
line width may be smaller than the 3D velocity dispersion, 
depending on the inclination and the amount of flattening. 
For example, for a flat BLR conflguration viewed at an in- 
clination angle i (where i = 90 corresponds to face-on), the 
line-of-sight velocity is smaller than the 3-D velocity by a 
factor of sin(«). The inferred mass would thus be smaller 
than the actual mass by a factor of sin^(i). For a ran- 
dom distribution of inclination angles, the average inferred 
mass would be decreased by this factor weighted by the dis- 
tribution, < sin^(i) >= An /J^^^ sin^(i) cos(i)di/47r = 1/3. 



Actually, the distribution is probably not random. In the 
unified scheme Seyfert 1 nuclei are viewed more face-on, 
within an opening angle of < 60 degrees, so that the 
weighted average mass-reduction factor is < sin^(i) >= 

47r J^^^sm^{i)cos{i)di/2n = 1/12. If the BLR distribu- 
tion is not flat but with an angular distribution of orbits 
spanning an angle S, the eff'ect would be less - between 
sin(z) and sm{d + i), depending on the distribution of or- 
bits, so that roughly < sin^(i + 5) > > 0.1. 

This scenario can explain the lower MBn/ibuigc values 
of NLANs. The distribution of inclination angles can also 
explain the larger scatter of NLANs, compared with broad 
line Seyferts. In this scenario, NLANs are AGNs with a 
flattened BLR geometry, seen nearly face-on. 

6.5. Uncertainties and non-virial BLR dynamics 

Is it possible that the black hole masses of NLANs are 
systematically underestimated by the reverberation - virial 
method? The reverberation masses of all AGN seem to 
be more or less correct, in light of the good agreement be- 
tween AGNs and quiescent galaxies in the case of the black 
hole - velocity dispersion relation. The uncertainty in the 
reverberation virial method is not well known. While for 
individual objects a factor of 2-3 may be representative, 
the sample average error is probably much smaller. This 
comes on top of the measurement uncertainty; the im- 
certainty in the BLR size (calculated from the cross cor- 
relation of the continuum and line light curves) can be 
quite large for objects with poor sampling or low variabil- 
ity. Note that if the virial assumption is incorrect and the 
gas in the line emitting gas in NLANs is imboimd, the 
gas velocity is actually larger than Keplerian, the black 
hole mass would have been overestimated, which would 
increase the BH-bulge discrepancy between NLANs and 
broad line Seyfert galaxies (and inactive galaxies). Simi- 
larly, if the BH mass of broad line AGNs were underesti- 
mated, their present good agreement with the BH-bulge 
relation of quiescent galaxies would deteriorate, as demon- 
strated in section 3.4. 

More generally, we may ask how reliable are the BH 
masses derived by the reverberation virial method? Could 
there be a systematic effect? The best resolution of HST 
in relatively nearby galaxies translates into a distance of a 
few tens of parsecs from the MBH. Reverberation mapping 
of the broad emission line region in AGNs gives a much 
closer view - a few light days from the center. Assuming 
the line-emitting matter is gravitationally bound, having 
a Keplerian velocity dispersion , it is possible to estimate 
the virial mass: M w G^^rv^. This expression may be 
approximately valid also in the case the line emitting gas 
is not bound, such as radiation-driven motions and disk- 
wind models (e.g. Murray et al. 1998). The main problem 
in estimating the virial mass from the emission-line data 
is to obtain a reliable estimate of the size of the BLR, 
and to correctly identify the line width with the velocity 
dispersion in the gas (WPM, Krolik 2001). WPM used 
the continuum/emission-line cross-correlation function to 
measure the responsivity-weighted radius ct of the BLR 
(Koratkar & Gaskell 1991), and the variable (rms) com- 
ponent of the spectrum to measure the velocity dispersion 
in the same part of the gas, which is used to calculate the 
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BLR size, automatically excluding constant features such 
as narrow emission lines and galactic absorption. Kaspi et 
al. (2000) find that for most objects in their sample, simi- 
lar mass estimates are obtained also using either the rms 
as or the mean FWHM. The case for a central MBH dom- 
inating the kinematics in the broad emission-line region is 
supported by the Keplerian velocity profile (z; a r~^' ^) de- 
tected in NGC 5548 (Peterson & Wandel 1999) and NCG 
7469 and 3C390.3 (Peterson & Wandel 2000). Although 
different emission- lines have quite different line widths, the 
delays also vary from line to line in such a way that the 
virial masses derived from different emission lines are all 
consistent with a single value. This demonstrates the case 
for a Keplerian velocity dispersion in the line- width/time- 
delay data. Another support for the validity of the rever- 
beration virial black hole mass estimates comes from the 
consistency of the black hole -bulge relations for AGN and 
inactive galaxies: the Mbh — c relation (Ferrarese et al. 
2001), and the Mbh — -^buigc relation (this work). 

6.6. Evolution 

The different BH-bulge relation of NLANs compared 
with broad-line AGNs and inactive galaxies could suggest 
that NLANs may be intrinsically different. One possibility 
is that they are in a different evolutionary stage. Wandel 
(1999) suggested that in AGNs which are radiating near 
the Eddington luminosity, the central black hole accretes 
within a relatively short time scale most of the readily 
available matter. BHs accreting near the Eddington limit 
are still growing, hence their mass (relative to the bulge) 
may be smaller than in AGNs with lower L/M ratios such 
as Seyferts Is which may be past the most active phase. 
Similarly, MBHs in normal galaxies may be the remnants 
of a past luminous AGN phase, and their BHs did not 
grow significantly after the decline in the accretion rate. 
This could explain the discrepancy between narrow and 
broad line AGNs (NLANs tend to be near the Edding- 
ton limit), as well as the similar BH-bulge relations of 
broad line AGNs and inactive galaxies. This scenario has 
been followed with a simple model calculation by Wang, 



Biermann and Wandel (2000) who examined the relative 
contribution of mergers to the bulge and BH mass. For 
reasonable parameters they find that a BH/bulge mass 
ratio of 0.002 is obtained by evolution of the MBH due 
to Eddington- limited accretion enhanced due to mergers. 

7. SUMMARY 

We re-examine the purported discrepancy between the 

black hole - bulge luminosity relation of MBHs in nearby 
galaxies (measured with stellar- and gas dynamics meth- 
ods) and reverberation -mapped Seyfert nuclei (Wandel 
1999). Using updated data, we find that the BH/bulge 
relation of AGNs and inactive galaxies are in good agree- 
ment over 3 orders of magnitude in bulge luminosity or 
BH mass. We show that the apparent difference reported 
previously originated from ovc;rc;stimated BH masses in 
quiescent galaxies, overestimated bulge luminosity for the 
Seyferts and including narrow-line Seyfert Is which seem 
to have genuinely low BH/bulge ratios. We find that the 
lower BH/bulge ratio of narrow line AGNs is part of a more 
general inverse correlation of the BH/bulge ratio with the 
emission-line width. Considering the BH-velocity disper- 
sion correlation, the lower BH/bulge ratio of narrow-line 
AGNs seems to be related (at least in part) to the host, as 
narrow and broad line AGNs (as well as inactive galaxies) 
seem to have a similar a— Mbh relation. We predict the ve- 
locity dispersion of high luminosity AGN using the Faber 
Jackson relation derived from Seyfert galaxies for with a 
measured cr. We discuss three classes of explanations for 
the lower BH/bulge ratios of narrow-line AGNs: black 
hole mass - bulge luminosity relation of Seyfert galaxies 
and MBHs in normal galaxies: Observational or method- 
related errors or bias, intrinsic and orientation-related ef- 
fects and Evolution. 
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Gebhardt, Luis Ho, Ari Laor and David Merritt. 



REFERENCES 



Bahcall, J.N., Kirhakos,S.,Saxe,D.H. & Schneider, D.P. 1997, ApJ, 

479, 642 Princeton U. 1987 
Cranshaw, D.M. 1986, ApJS62, 861. 
Erackleous, M. and Halpcrn, J. P. 1994 ApJS, 90, 1 
Ferrarese, L. and Merritt, D. 2000 ApJ, 539, L9 
Ferrarese, L., Pogge, R.W., Peterson, B.M. , Merritt, D. , Wandel, 

A. & Joseph, C.L. 2001, ApJ, 555, L79 
Gebhardt, K. et al. 2000a ApJ, 539 L13 
2000b ApJ, 543, L5 

Hamilton, T.S., Casertano, S. & Turnshek,D.A. 2001, ApJ, in press 

(astro-ph 0011255) 
Ho, L.C. 1999 in Observational Evidence for Black Holes in the 

Universe, ed. S.K.Chakrabatri (Dordrecht: Kluwer) 153 
Kaspi, S. et al. 2000, ApJ, 533, 631 

Kirhakos, S., Bahcall, J.N., Schneider, D.P. & Kristian, J. 1999, 
ApJ, 520, 67 

Koratkar, A.P., & Gaskell, CM. 1991, ApJS, 75, 719 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Kormendy, J. & Gebhardt, K. 2001, proc. of the 20th 
Texas Symposium, eds. H.Maxtel & J.C.Wheeler, in press 

(astro-ph/0105230) 
Krolik, J.H. 2001, ApJ, 551, 72. 
Laor, A. 1998, ApJ, 505, L83 
Laor, A. 2001, ApJ, 553, 677 
Magorrian, J. et al. 1998, AJ, 115, 2285 



Mathur, S.,Kuraszkiewicz, J. &; Czerney, B. 2001, New Astronomy, 

in press (astro-ph/0102356) 
McLure, R.J. & Dunlop, J.S 2001, MNRAS, in press (astro- 

ph/0009406, MDOl) 
Merritt, D. and Ferrarese, L. 2001a, MNRAS, 320, L30 

2001b ApJ, 547, 140 

Murray, N. 1998, ApJ, 494, 125 
Nelson, C.H. & Whittle, M. 1995 ApJS, 99 67 
Nelson, C.H. 2000, ApJ, 544, L91 
Neugebauer, G. et al. 1987, ApJS, 63, 615 
Peterson, B. & Wandel, A. 1999, ApJ, 521 L95 

2000 ApJ, 540, L13 

Richstone, D. et al. 1998, Nature, 395, A14 

Simien, F. & de Vaucouleurs, G 1986, ApJ, 302, 564 

Santos Lleo, M., Clavel J., Schulz, B. et al. 2001 Astronomy and 

Astrophysicsin press (astro-ph0102356) 
Virani, S.N., De Robertis, M.M. and Van Delfsen, M.L. 2000, 

AJ120 1739 
Wandel, A. 1999, ApJ, 519, L39 

2000, in "Probing the Physics of AGN by Multiwavclength 

Monitoring", ASP 224, 365, eds. B.M. Peterson.P.S.Polidan & 
R.W.Pogge 

2001, BAAS, 33, 898. 

Wandel, A. & Mushotzky, R.F. 1986, ApJ, 306, L61 
Wandel, A., Peterson, B. and Malkan, M 1999 (WPM), ApJ, 526, 
579 



12 



Wandel 



Whittle, M. et al. 1992, ApJS, 79, 49 



Wang, Y.P., Biermann, P. & Wandel, A. 2000, Astronomy and 
Astrophysics, 361, 550 



